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A Problem in Self-Heating of a Spherical Body 

S. M. Genensky* 

An analytic steady-state solution is developed for a spherical body in which heat is 
generated according to a first-order unimolecular-reaction law, and lost at the surface in 
accordance with Newton's law of cooling. The temperature within the sphere depends 
chiefly on the radial distance from the center, but also on the ambient temperature, the 
surface heat-loss coefficient, and the material properties, which are assumed constant. 



1. Introduction 

The analysis described was pursued in connection 
with an investigation into the self-heating of fibrous 
materials carried on at the National Bureau of 
Standards. 1 

Comparison of the results of this analysis with the 
experimental results indicated that the assumption of 
a single first-order unimolecular reaction was an 
oversimplification for the material investigated 
experimentally. However, the analysis serves to 
indicate the relative importance of the various 
parameters entering into the problem, and may be 
helpful in indicating the kind of data that would be 
of greatest use in further investigations in which the 
assumptions may be applicable. 

A homogeneous sphere of radius B is generating 
heat, under steady-state conditions, in accordance 
with a first-order unimolecular-reaction law, Ae~ {E/RT \ 
and losing heat from its surface to an atmosphere at 
constant temperature, T a , according to Newton's 
law of cooling. The temperature of the sphere, T(r), 
is independent of time and depends upon the radial 
distance, r, from the center of the sphere. Within 
the sphere, transfer of heat by means other than 
conduction is considered negligible. The thermal 
conductivity, k; gas constant, R; rate of reaction- 
frequency product, A; and activation energy, E, 
of the material of the sphere are assumed to be known 
constants, as are the temperature, T , and gradient 
dT/dr=0, at the center of the sphere. The problem 
is to find the surface temperature of the sphere, 
T B , at r=B, and the heat-transfer coefficient, h, 
between the surface of the sphere and a surrounding 
atmosphere. 

2. Analysis 

Mathematically the problem becomes 
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Then eq (1) to (4) become 
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It is further assumed that (a) the temperature 
gradient within the sphere exists and is continuous 
throughout the interval 0<t;<1, and (b) that (he 
absolute value of the difference between T(rj), 
0<77<1, and T is small in comparison to T . Under 
these restrictions, a technique presented by Cham- 
bre 2 proves useful to observe that 



V Vo-AV^Voi 1 ^ 1 vJJ-vSAVo)' 



(9) 



where AV=V —V. Now, because |AF/F |<1, the 
series converges, and for a sufficiently large positive 

iV 

integer, N, (l/V )J2(AV/V ) i is a very good approx- 

imation to 1/V. 

Therefore, (5) may be written approximately as 



N 



dr) 2± r)dri^ 



=0 (0<)<1). (10) 



2 P. L. Chambre, On the solution of the Poisson-Boltzman equation with ap- 
plication to the theory of thermal explosions, J. Chem. Phys. 20, 1795 (1952). 
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and in the special case V= V 

{N+l)-B l V+B 2 V 2 -. . .+(-l) N B N V N =l 



N=0,l,2, . . 



Therefore, (10) may be written 






(ii) 



Because (11) is analytic and regular except at r)=0, 
and further as dV/dr)=0 at 97=0, the solution of (11) 
may be written in the form 

V=V(r ] )=^Za i v i (0 < v < 1). (12) 

i=0 

Recalling (6), eq (12) at rj=0 yields 

V =V(0)=?jlr==a , (13) 

and further differentiating (12) once and considering 
(7) , it is found that 

a x = 0. (14) 

Differentiating (12) twice, substituting these first 
two derivatives into eq (11) and rearranging terms: 
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where D=Ce Vo and the a t are given by 
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,(15) 
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ai ~ (i 2 +i){i-2)\' 
for 17=0, and i=2,3,4, 



drj'- 



(16) 



Using the coefficients obtained from (13), (14), 
and (16) and substituting them into eq (12), V can 
be computed for 0<ij<1. 

In particular, for 77= 1, eq (12) gives the dimension- 
less surface temperature V(l), and thus recalling 
that T(rj)=(EIR)V(ri) } the surface temperature, T(l), 
which equals T B , is easily found. Differentiating 
eq (12) once and evaluating this derivative at 77=1, 
h can then be found by solving eq (8), for all the 
other factors involved in this equation are either 
given or have been computed. 



The first six nonzero coefficients have been evalu- 
ated and are 
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Observe that a 2 *-i (i= 1,2,3, . . .) are zero. This 
follows from ai = and the spherical symmetry of 
the problem. 

Figure 1 is a plot of V(l), as a function of C and 
1/Vo, where C was allowed to vary over the range 
10 -1 <C<10 19 and 1/V took on selected values in 
the range 5<l/Vo<100. The points were obtained 
by using 

10 
V(v)=^2a iV i (17) 



over its range of applicability for the values of C 
and 1/Vo within the above limits. Ipr^i 

In the foregoing, it has been assumed that the 
temperature at the center of the sphere T is known, 
and also that the physical condition indicated by 
eq (8) is satisfied. 

The more common steady-state problem of finding 
the temperature at the center of a sphere that gener- 
ates heat according to a first-order unimolecular 
reaction law and loses heat according to Newton's 
law of cooling may be solved by using the analysis 
described above. In this case, E, R, k, A, h, T a , 
T B , and B are assumed known. 

If the sphere being considered really satisfies the 
assumptions of the problem, then its temperature 
mast satisfy eq (1). Because eq (1) is of the second 
order, only two boundary conditions can be pre- 
assigned. Further, because the temperature at the 
center of the sphere must remain finite, one of these 
preassigned boundary conditions must be dT/dr=0 
at r=0. This leaves but one free boundary con- 
dition. However, at the surface of the sphere both 



(f)„,— t {T -- Ti < is > 



and 



T=T B at r=B (19) 

must be satisfied ; and because only one of these two 
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Figure 1. V (1) as a function of C for selected values of 1/V Q . 



boundary conditions may be formally assigned, it 
must be assumed that the other is consistent with 
the problem. Therefore, figure 1, which represents 
eq (17), may be used to solve this problem. Because 
R, E, k, A, and B are known, both C and V(l) (recall 
that T B =EV(1)/R) can be computed. Thus, by 
using figure 1, l/V may be found by interpolation, 
provided the point (C, V{\)) under consideration 
lies within the region of applicability of eq (17), as 
indicated above. Because the center temperature 
of the sphere is T , that is EV /R, the problem is 
solved. 

Here no use was made of condition (18), and as 
mentioned earlier, this condition must be satisfied 
independently, if the assumptions made are fulfilled. 
Thus the degree of consistency between conditions 
(18) and (19) serves as an indication of the applica- 
bility of the assumptions made in this analysis. 



Washington, March 19, 1954. 
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